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NOMENCLATURE 


a Width of the DHE duct, m, 

b Height of DHE duct, m, 

A Area, m 2 . 

Ahp Cross section flow area of vapor in heat pipe, m 2 

C Specific heat, J/kg K, 

D Fuel pin, or pipe diameter, m, 

D e Equivalent hydraulic diameter, m, 

f Friction coefficient, 

F Radiation view factor, 

g Acceleration of gravity, m/sec 2 , 

h Heat transfer coefficient, W/m 2 K, latent heat, J/kg, 

H Height of reactor, m 

hfg Latent heat of vaporization, kJ/kg, 

k Thermal conductivity, W/m K, 

L Length of the coolant flow path in DHRL, m 
m Coolant mass flow rate, kg/s, 

Nu Nusselt number, (h De/k ), 

Nu s Slug Nusselt number, equation (24), 
p Pressure, MPa, 

P Pressure, MPa, or power, W, or pitch distance, m, 
P e Peclet number (Re. Pr), 

Pr Prandtl number, (C p p/k), 

P/D Pitch-to-diameter ratio, 

Volumetric decay heat , W/m 3 , 
qi Linear power, W/m, 

r Radial coordinate, m, 

R Radius of fuel pin, m, 

Re Reynolds number , (pvDe/p), 
s Diameter of the wrap wire, m, 

t Time, s, 

T Temperature , K, 


V 



U Overall heat transfer coefficient, Equation. (22), 

v Flow velocity, m/s, 

z Distance from the entrance of reactor core, m, 

Z Dimensionless distance from the entrance of reactor core, m. 

Subscripts 

a Ambient, 

b Coolant bulk, 

c Lithium coolant, cladding, 

eff Effective 

f Flow, fuel, 

g Fuel-cladding gap, 

h Convection heat transfer, heat transfer, 

L Liquid coolant, DHRL, 

hp Heat pipe, 

i Inner surface, 

m Log-mean, lunar surface, 

o Initial, 

r Radiation, 

rej Rejection, 

loss Losses, 

sonic At heat pipe sonic limit, 

s Cross-section, 

v Saturated vapor, 

w Wall. 

Superscripts 
Core SP-100 reactor core 

Pipe DHRL pipes 

Greek Letters 
8 Thickness, m, 

e Emissivity, 

a Stefan-Boltzmann constant 

(5.669x10-8 w/m 2 K 4 ) 
p Dynamic viscosity, kg/m s, 

p Density, kg/m 3 . 
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1. INTRODUCTION 


A permanent outpost on the moon is being considered by NASA for the 
early part of the next century. The lunar reference mission objectives are 
prioritized as follows: (1) to establish a lunar outpost with a long-term manned 
presence, (2) to gain experience in working on planetary surfaces, (3) to conduct 
scientific and manufacturing resources, and (4) to develop crew planning 
capabilities. The operation concept for the lunar outposts provide for the 
eventual decentralization of operation control, whereby planning, monitoring and 
controlling functions will be shifted from being earth-based to being shared by 
earth and the outpost. T 2 The power requirements and decentralization of 
operation control will evolve within mission’s successive phases. These phases are 
precursor/emplacement, consolidation and utilization/demonstration. While the 
power requirements for the precursor/emplacement phase are in the tens of 
kilowatts electric, they could grow to the hundreds of kilowatts electric in the 
utilization/demonstration phase. 

Operation analysis of the manpower, machine resources and power 
requirements through the year 2,013 for the lunar program was conducted.! ,2 
The stationary power systems for the lunar outpost will be designed to meet the 
evolutionary growth in power demand ranging from the tens to hundreds of 
kilowatts electric. As the power demand increases for a constructible habitat, an 
SP-100 reactor, coupled to a number of dynamic energy conversion engines, is 
emplaced to supply 550 kW e . This power level has been selected based on the 

NASA Space Exploration Initiative 90-day study. 1 

A recent study suggested that an SP-100, 550-kWe power system, with four 
Stirling or Brayton engines (three operating and one standby) would be optimum 
for reliability consideration. 4 In addition to system reliability and redundancy, it 
is important that the integrity of the nuclear reactor core be maintained during a 
loss-of-flow accident. As a result, the power system is designed with passive 
means for the removal of decay heat by natural circulation from the reactor core 
after shutdown. The decay heat removal loop of the system is equipped with a 
separate heat exchanger and a heat rejection radiator on the lunar surface. More 
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details on the power system description and operation are included in the next 
section. 


This research investigated the decay heat removal from the SP-100 reactor 
core of a 550-kWe power system for a lunar outpost by natural circulation of 
lithium coolant. A transient model that simulates the decay heat removal loop 
(DHRL) of the power system was developed and used to assess the system s decay 
heat removal capability. The effects of the surface area of the decay heat rejection 
radiator, the dimensions of the decay heat exchanger (DHE) flow duct, the 
elevation of the DHE, and the diameter of the rise and down pipes in the DHRL 
on the decay heat removal capability were examined. Also, to determine the 
applicability of test results at earth gravity to actual system performance on the 
lunar surface, the effect of the gravity constant (lg and l/6g) on the thermal 
behavior of the system after shutdown was investigated. 

The following section presents background information on the design and 
the materials used for the SP-100 nuclear reactor and a detailed description of the 
SP-100 fuel element. Section 3 gives a detailed description and discussion of the 
NASA SP-100, 550-kWe power system for a lunar outpost. Section 4 presents 
the mathematical formulations, governing equations, and method of solutions for 
each of the component's model as well as for the fully integrated decay heat 
removal loop model. Section 5 presents the transient results on the decay heat 
removal capability of the DHRL with emphasis on both time and spatial changes 
of the coolant temperature within the DHRL. In this section, the results of the 
parametric analyses, investigating the effects of DHE dimension and elevation, 
decay heat radiator area, guard vessel heat pipe diameter, external pipe diameter 
and gravity on the passive removal of decay heat in the SP-100, 550-kWe power 
system, are also presented and discussed. In section 6, a summary of the 
important results is given and conclusions are stated. More information on the 
mathematical formulations is given in the appendices at the end of the report. 
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2. BACKGROUND 


The SP-100 space nuclear power system technology is being developed 
under a joint program between the Department of Energy (DOE), National 
Aeronautics and Space Administration (NASA), and the Department of Defense 
(DOD) to support future space missions with electric power requirements in the 
tens to hundreds of kilowatt electric.3 The base 100-kWe power system design 
for a space mission to be developed in the earth s orbit utilizes SiGe 
thermoelectric (TE) generators to partially convert the reactor thermal power 
into electrical power with an efficiency of about 6 percent. To enhance the 
electrical power output of the SP-100 system to supply surface power for a lunar 
base, NASA will replace the TE converter with a dynamic system, such as free- 
piston Stirling (FPS) or Closed Brayton Cycle (CBC) engines, with an efficiency 
that could be more than three times larger. Therefore the SP-100 nuclear reactor 
technology and design can be used, while the electrical output can be increased to 
more than 500 kWe by such a power system. 

The SP-100 nuclear assembly includes the reactor, shield and control 
systems. The SP-100 is a fast spectrum reactor cooled by pumped liquid lithium. 
The critical components of the reactor are: (1) fuel and cladding, (2) honeycomb 
structure, (3) reflector supports, (4) safety rod thimbles, (5) grid for pin support, 
(6) reactor vessel, (7) core support structure, and, (8) piping that interfaces with 
the primary heat transport subsystem. Figure (1) displays these components in 
the flight configuration and the core layout. The reactor core, in a triangular 
arrangement, is enclosed in a niobium refractory alloy known as a PWC-11 
vessel. This alloy is niobium - 1% zirconium with 1000 ppm of carbon added to 
improve its high temperature creep strength. 

The fuel pins are arranged in the reactor in approximately 61 hexagon 
assembly groups. Each group is housed with one of the hexagonal cells with the 
honeycomb structure of the core. This honeycomb structure prevents 
unacceptable fuel pin blowing and fluid cross flowing within the reactor core. A 
shadow shield (LiFl/W) that attenuates both neutrons and gammas is employed to 
protect the spacecraft electronics from the radiation emitted for the reactor. The 
same reactor design and shadow shield assembly will be used in the 550-kWg 
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power system for a lunar outpost. Given a reactor thermal power of 
approximately 2.3 MWt and nominal conversion efficiency of the dynamic 
conversion engines of 23 %- 25 %, a net electrical power output of 529-575 kWe 
can be obtained. 

The fuel pin structure is detailed in Fig. (2). Because uranium nitrite (UN) 
has high thermal conductivity, high uranium density and low fission gas 
release/swelling, the fuel pins in the SP-100 utilize it as fuel material. The UN 
fuel region in the fuel pin is bounded by Berylium Oxide neutron reflector 
materials on both ends. A plenum is located at the top of the upper reflector to 
accommodate the fission gasses released from the fuel. A rhenium (Re) liner is 
placed between the helium filled fuel gap and the Nb-1% Zr cladding. This Re 
liner is used to enhance flooding subcriticality and high temperature strength and 
to provide a chemical barrier. On the outer surface of the cladding, a wire wrap 
is attached to maintain proper space among fuel pins and to enhance heat transfer 
from the fuel pin to the coolant during normal operation. The fuel pins are 
arranged in a triangular lattice with a pitch-to-diameter ratio of 1 .07 . 
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Fig. 2 Fuel Pin Configuration. 



3. SP-100, 550-kWe POWER SYSTEM FOR A LUNAR OUTPOST 

The SP-100, 550-kWe nuclear power system for a lunar outpost is 
currently being developed by Rockwell International Corporation for the 
National Aeronautical and Space Administration (NASA). 4 The power system 
employs an SP-100 nuclear reactor cooled with liquid lithium, which is coupled 
to four dynamic energy conversion engines (either a Brayton or Stirling power 
conversion subsystem). Figure (3) shows a pictorial view of the SP-100, 550 
kWe power system for a lunar outpost. The whole system is emplaced inside a 
cylindrical, excavated cavity in the lunar regolith. A cavity liner (Guard Vessel) 
is provided as an integral part of the system. In the primary loop, the lithium 
coolant from the core is circulated by an electromagnetic (EM) pump, through 
the primary side of an intermediate heat exchanger, then returned to the reactor 
core. An expansion tank with a free lithium surface is located at the high point of 
the primary loop. This allows for expansion of the lithium as it heats up and for 
collection of the helium gas formed by the neutron interaction with the lithium. 
The primary loop is thermally coupled, via the intermediate heat exchanger, to 
either liquid lithium (for Stirling System ) or He-Xe gas (for Brayton System) in 
four secondary loops. Each is equipped with a pump and a dynamic energy 
conversion (Stirling/Brayton) engine. The intermediate heat exchanger isolates 
the power conversion heat transfer fluid from the reactor coolant, thereby 
reducing the possibility of fission product contamination of any component 
located at the lunar surface. The primary loop EM pump is a flat linear induction 
pump with redundant starters and power supplies. 4 

As shown in Fig. (3), the guard vessel fits closely around the reactor and 
the primary coolant loop to ensure that the reactor core is always covered with 
lithium coolant in the event of a small leakage happening in the primary loop or a 
partial/total loss of flow. The entire power system, with the exception of the 
radiators, is emplaced in a cylindrical, excavated cavity in the lunar regolith. 
The primary coolant loop, the radiation shadow shield and the SP-100 reactor are 
located at the bottom of the cavity and the secondary loop components are placed 
at the top of the cavity. The cavity is lined with a stainless-steel guard vessel, 
which is an integral part of the system and is cooled on its outer surface by heat 
pipes. The guard vessel heat pipes are integrated with those of the dedicated 
radiators on the lunar surface for auxiliary and decay heat removal. 
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These decay heat radiators are located between the main radiator panels. In 
addition to the obvious structural and integrational advantages, the guard vessel is 
used to transport the heat losses from the reactor/primary cooling system and the 
energy deposited in the vessel structure by neutrons and gammas during nominal 
power operation, as well as the reactor decay heat after shutdown. 

The flow schematic and state points for the Brayton system are shown in 
Fig. (4). A single primary Li loop transports heat from the reactor to the 
Brayton primary heat exchanger. The Brayton cycle uses an He-Xe gas mixture 
as the working fluid. Heat is rejected from the cycle through a NaK heat 
rejection loop. Each Brayton loop is cross-coupled to each of four radiator 
panels so that if a power conversion failure occurs, there is no loss of radiator 
area. The NaK loop is also used to cool the Brayton alternator. 

The flow schematic and state points for the Stirling system are shown in 
Fig. (5). The primary and secondary Li loops are used to transport heat to the 
Stirling engine heater. Helium is used as the working fluid for the Stirling cycle. 
Heat is rejected from the Stirling cooler through a Nak heat rejection loop, which 
is also used to cool the alternator. As with the Brayton system, the heat rejection 
loops in the Stirling system are cross-coupled to each radiator panel. 

Figure (6) presents a line diagram of the overall deployed arrangement for 
the Stirling/Brayton system with the exception the guard vessel and decay heat 
radiators. As shown in Fig. (6), the four conversion engines are coupled in 
parallel to four separate radiator panels on the lunar surface. Therefore, a failure 
of one of the engines or a loss-of-flow in one of the secondary loops will not 
affect the operation of the other three engines. In this case, the heat rejection of 
the system will be handled by the four radiator panels, but only three engines 
would supply the 550 kWe needed for the lunar outpost operation. 

As Figs. (3) and (6) show, the decay heat exchanger (DHE) is an integral 
part of the primary coolant loop. The DHE is hydrodynamically coupled to the 
primary coolant loop in the hot line and in the return line via a venturi. The 
DHE has a rectangular flow duct and, as shown in Fig. (3), it is equipped with 
small sodium heat pipes to increase its effective surface area. There are about 100 
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Fig. 5 SP- 100,550 kWe Stirling Flow Schematic and State Point. 
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such sodium heat pipes measuring 30 cm in length and 2.75 cm in diameter in the 

current design. 22 

Figure (7) presents a schematic diagram of the DHRL. During normal 
power operation, the venturi in the primary coolant loop equalizes the pressure 
drop across the DHE hence net flow occurs through the decay heat exchanger is 
zero. However, following a hypothetical Loss-of-Flow Accident (LOFA), force 
flow through the venturi drops to zero, hence allowing natural circulation of 
lithium coolant through the DHE duct under the effect of the difference m 
coolant densities in the rise and down pipes of the DHRL (see Fig. (7». e eat 
rejection from the DHE is accomplished by radiation to the inside wall of the 
guard vessel, where heat is transported by the guard vessel heat pipes to the decay 
heat radiators on the lunar surface. These radiators are separated from those used 
for heat rejection from the energy conversion engines during normal full power 
operation of the power system. All primary loop components are electrically 
trace heated to provide for controlled Li thaw during start-up. 

To investigate the thermal behavior of the SP-100, 550-kWe power system 
for a lunar outpost after reactor shutdown, an integrated model of the system s 
decay heat removal loop by natural circulation of primary Lithium coolant has 
been developed. This integrated model is described in detail in the following 

section. 
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A Schematic Diagram of Decay Heat Removal Loop. 




4. MODEL DESCRIPTION 


The integrated decay heat removal model of the SP-100, 550- kWe system 
developed herein consists of three coupled sub-models: (a) two-dimensional 
transient thermal model of the fuel pin, (b) transient natural circulation thermal- 
hydraulic model of the DHRL, and (c) heat rejection model for the DHE/guard 
vessel heat pipes. The properties of the lithium coolant and those of the fuel, 
cladding, and structural materials are taken to be temperature dependent. As 
indicated earlier, the decay heat exchanger is radiatively coupled to the inside 
surface of the guard vessel. The radiation view factor between the DHE and the 
guard vessel is calculated as a function of the dimensions of both the DHE and the 
guard vessel. 6 (In Appendix A, an analysis is performed to estimate the radiation 
view factor between the DHE and guard vessel.) However, the effective surface 
emissivity for the DHE wall and the inside surface of the guard vessel are 
assumed constant during the analysis (see Table 1). Fig. (8) shows a block 
diagram of the DHRL model. 

4.1 Two-Dimensional Transient Thermal Model of Fuel Pin 

The fuel pin model calculates the radial and axial temperature distributions 
inside the fuel pin as functions of heat generation in the fuel (fission power or 
decay power), coolant mass flow rate, and coolant temperature in the reactor 
core. Although the two-dimensional heat transfer equation can be used for this 
purpose, this can be simplified by neglecting the axial heat conduction. This 
assumption is justified by the negligible axial heat conduction due to the long 
length of the fuel pin compared to its radius and high thermal conductivity of the 
UN fuel. In this model, the fuel pin is discretized into small axial segments (see 
Fig. (9)), and in each segment the radial transient heat conduction equation is 
solved while the axial heat conduction between segments is neglected. The axial 
temperature distribution in the fuel pin is obtained by thermal-hydraulically 
coupling the different axial segments. 
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4.1.1 Governing Equations 

In each axial segment in the fuel pin, only the radial heat conduction is 
considered. The transient radial heat conduction equation is: 


pf Q 


W=il(k f r 

at r ar v 


3T(r,t) 

ar 


) + q (r,t) 


( 1 ) 


In the cladding and the fuel-cladding gap, the second term on the right hand side 
of Eq. (1), the volumetric heat source, is zero. 

The heat source term q'"(r,t) in Eq. (1) is the decay heat inside the fuel 
region in our decay heat removal calculation. This decay power is the function of 
the fission power level and reactor operation time before shutdown and fuel 
material composition. In our analyses, the decay heat generation in the fuel 
region is calculated using the decay heat curves recommended by Marr and 
Bunch. 10 Fig. (10) presents these decay heat curves for fast spectrum reactors as 
a function of time after shutdown. These curves are plotted using a long reactor 
operation time before shutdown (> 6 months). 

The initial and boundary conditions for the fuel pin model are: 

T(r,0) = T o (r),(«| = 0 (2a) 

k, (^r i ) =h(mt) ' Tb(t)| (2w 

In Eq. (2b), the heat transfer coefficient between the lithium coolant and 
the cladding, h, is calculated using a natural convection heat transfer correlation 
for triangular lattices bundles:^ 

Nu = 24.0(( Af/A h ) 18 . Pe)° 6 , (3) 
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Fig. 10 Decay Heat Curves for Fast Spectrum Reactors. 


where Af and Ah are the flow area and the heat transfer area in the core 
respectively. 

In the gap between fuel and cladding, the gap thermal conductance is 
calculated using either a closed or an open gap model, which allows for fission 
gas release into the gap during reactor operation. 7 For the open gap case, the 
gap conductance is given as: 


h = h k + h r . 


(4) 


In this equation, h k is the contribution of thermal conduction through the 
gas mixture inside the gap and h r is the thermal radiation contribution of between 
the surface of the fuel and the cladding. 

For the closed gap, the overall gap conductance is taken as the sum of 
conduction and radiation contribution through the gas entrapped between the fuel 
and cladding inner surface, hk and h r , and the solid-solid conduction at the contact 
point, h s : 


h — h k + hf + h s . 


(5) 


The thermal conduction through the gas gap is given as: 18 


h t = 


(5f + S c ) + G + (gf + g c ) 


( 6 ) 


where S is the surface roughness, g is the jump distance and G is the gap 
thickness. The subscripts f and c represent the fuel and cladding respectively. In 
Eq. (6), the thermal conductivity of the gas mixture, k g , is calculated as: 


k g = (k,p(k2p(k 3 ) Z3 


(7a) 


where 


ki = Aj T g ° 79 , (i=l ,2,3) 


(7b) 
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In Eq. (7b), T g is the average temperature of the gap; A i 1S 15.8x10 6 , 
0.72x1 0 6 and 1.15xl0- 6 for He, and for fission gases Xe and Kr respectively, Xi 
is the mole fraction of a particular gas in the mixture. In Eq. (6), the gap 
thickness, G, becomes zero when this equation is used in the closed gap equation. 
For the closed gap case, the thermal conductance due to solid to solid contact is 

given as: 



where C is a constant (0.01), P is the interfacial pressure (MP a ) and H is the 
hardness of the softest material, which is Nb-l%Zr cladding. The effective gap 
thickness for closed gap, 5 e ff, is given as: 


8etf=Y(^^> (9) 

Finally, the radiation contribution is given as: 

h, = ae e ( T f 2 + T c 2 ) (T f + T c ), (l0) 

where a is the Stefan-Boltzman constant (5.67x10-* W/m 2 K 4 ), the effective 
surface emissivity, £e is calculated by: 


e e = 


1 r f i 

-+—(—- 1 ) 

£f r. e r 


-1 


(ID 


The volumetric heat generation in the fuel pin is allowed to vary in a 
cosine distribution in the axial direction. According to El-Wakil, 19 the volumetric 
heat generation at the geometrical center of a heterogeneous core, q C o can be 
calculated by: 
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( 12 ) 


Qco — 


Pf/1.05 
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where Pp, is the total fission power (or decay power) generated in the fuel 
region, n is the number of the fuel pins in the reactor, A s is the cross section area 
of the fuel pellet in the pins, and He and Re represent the extrapolated height and 
radius of the core. Since the normal neutron flux distribution is assumed in Eq. 
(12), the hottest fuel pin in the core is located at the center of the cylindrical core 
at r=0. The axial distribution of the volumetric heat generation for the fuel pin in 
the hottest channel at the center of the reactor core can be given as: 


q ,n (z) = qcocos(2£) 
We 


(13) 


The mass flow rate and coolant temperature in the reactor core are 
determined from the coupling of the fuel pin thermal model with the thermal- 
hydraulic model of the DHRL. 


4.1.2 Method of Solution 


The transient heat conduction equation (Eq. (1)), together with the 
boundary conditions (Eqs. (2a) and (2b)) and all the input parameters (Eqs. (3)- 
(13)), is solved in each axial segment of the fuel pin using a finite element method 
combined with an implicit time integrator.^ 


In each segment of a fuel pin, the transient one-dimensional heat 
conduction equation can be written as: 


Pf Cf 


w.iite?)-,- 

dt r 3r 3r 


= 0 


(14a) 


Defining: 


[T(r,t)]= Pt Cff-^(k fr f )-q' 


(14b) 
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and using Gerlerkin Approximation, 21 Eq. (14a) can be written as: 



NiL[ T*(r,t) ] dr = 0 

* 


( 15 ) 


where T*(r,t) is the approximation of T(r,t) is the solution domain of [0,R] which 
is defined as: 


M 

T*(r,t) = £ Nj(r) Tj(t) 


i=i 


(16) 


In Eq. (16), Nj is the basis function. 21 Substituting Eqs. (14a), (14b) and 
(16) into Eq. (15) and rearranging the result, yields: 



Equation (17) is equivalent to Eq. (14) under the Gerlerkin 
Approximation. 21 By choosing the basis function Ni(r), we can solve Eq. (17) 
for Tj(t), and then the approximate solution of T(r,t) in [0,R] domain, T*(r,t), 
can be calculated by Eq. (16). In Eq. (17), M is the number of elements used in 
the [0,R] domain. Finally, the matrix form of Eq. (17) can be written as: 

[C] {T} = [K] {T} +{Q}, (18) 
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where 


Ci,i Ci, 2 •••• c i.m 

C2.1 C2,2 ^2,M 

Cm,i Cm, 2 •••• Cm.m 


Ku K 1)2 
K2,l K2,2 

Km,1 Km,2 


1 Qi 
Q2 

( q ) = | : 

Qm-i 
, Qm 




and the coefficients in these matrices are defined as follows: 
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In Eq. (18), [C] is the capacitance matrix, [K] is the conductivity matrix, 
and {Q } is the heat source vector. The time integrator of the transient term on 
the right-hand side of Eq. (18) can be expressed as: 

(T} n+1 = {T}"+ s[a{f }"*' + (l-a)tf )"] , (19) 

where a is a constant determining the integration method and s is the time step 
size. As we can see in Eq. (19), a (which varies from zero-unity) represents the 
weight used for (T} n+ l and {T} n to calculate {T} n+ T The integration method is 
called 'fully explicit' for a = 0 and 'fully implicit' for a = 1.0. 


Using Eq. (18) in both nth and (n+l)th time steps, the following final form 
of the matrix equation is obtained as: 

( [C] - sa [K] ) {T}"+i = ( [C] + s(l-a)[K] ) {T}"+ s{Q}"+i , (20) 

Solving this equation with the appropriate matrix decomposing algorithm gives 
the time dependent, two-dimensional temperature distributions inside the fuel pin 
for the given core coolant temperatures, coolant mass flow rate, and fission (or 
decay) power. 
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4.1.3 Model Verification 


To verify the 2-D fuel pin, the calculated linear power distribution at 
steady-state is compared with the analytical value using in the equation. 

J rTm 

kf(T) dT 

Ts . (2D 

where T m and T s are the calculated temperatures at the fuel centerline and the 
cladding surface respectively. For a cosine power profile P(z) = P 0 cos(7tz/H e ), 
we calculated the linear power of the fuel pin using qi = P(z)/H e . These values are 
compared with the results calculated from Eq.(21) for a different mesh number. 
T m and T s are calculated using the 2-D fuel pin model at the same cosine power 
profile for UN fuel material. This comparison is presented in Figs. (11) and (12) 
for (10x10 = 100) meshes and (40x40 = 1600) meshes respectively. As these two 
figures show, for the 100 meshes case, the maximum difference between the 
numerical and analytical solution was 2.7% and for 1600 meshes was 0.044%. 
Increasing the mesh number increases the accuracy of the numerical results. 
However, 100 meshes still gives a reasonable accuracy and need much less CPU 
time for calculation compared with 1600 meshes. Therefore, in the subsequent 
analyses, a total of 10x10 meshes were used in the radial and axial direction of 
the fuel pin, respectively, to calculate its corresponding transient temperature 
profile. 

4.2 Thermal-Hydraulics Model of the DHRL 

The natural circulation thermal-hydraulic model of the DHRL (Figs. (7) 
and (8)) couples the reactor core, rise pipe (adiabatic section), decay heat 
exchanger (DHE), and down pipe and reactor downcomer (adiabatic sections). 
The energy and momentum balance equations are formulated and solved for the 
coolant bulk temperature within each region. These temperatures are used to 
update the coolant properties in the DHRL then, the overall momentum balance 
equation is solved for the coolant mass flow rate. This process continued, 
marching in time to calculate the coolant temperature and mass flow rate as a 
function of time after reactor shutdown. 
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4.2.1 Thermal Model 


From the energy balance in a control volume of the DHRL, the general 
spatial distribution of the coolant bulk temperature in the DHRL can be 
determined by: 



( 22 ) 


On the right hand side of Eq. (22), the first term and the third term 
correspond to axial convection and conduction heat transfer, respectively, and the 
second term is the radial heat transfer through the wall. In the reactor core, DHE 
and the adiabatic sections (rise/down pipes, reactor downcomer), the overall heat 
transfer coefficient, U, is calculated by different equations. In the core coolant 
flow channel, U is equal to hi, which is calculated by Eq. (3), and in the DHE 
duct, U is determined as: 

U = {X [ (hj Aj)’ 1 + (S^AJcJ -Kh, Aj.)' 1 ]} \ (23) 

Where, the convective heat transfer coefficient, hi is determined from the 
relations: 11 

Nu = 4.367, for Re < 2,100 , (24a) 

and, 

Nu =(2/3)Nu s +0.015 P e ° 8 , for Re > 3,000 . (24b) 

For 2,100 ^ Re < 3,000, Nu is determined from a linear interpretation of 
Eqs. (24a) and (24b). The slug Nusselt number in Eq. (24b), Nu s decreases with 
the aspect ratio of the DHE flow duct (a/b) increase (see Fig. (13)). 11 For the 
base case aspect ratio of 0.26 (see Table 1.), Nu s = 6.5. 11 
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The radiative heat transfer coefficient between the outer surface of the 
DHE wall and the inside surface of the guard vessel is given as: 

hf = eoF (Tw + Th P ) (T w + T hp ) . (25) 

The temperature of the guard vessel wall is assumed to be equal to that of 
the evaporator of the vessel’s heat pipes, T hp . The reactor decay heat is removed 
by thermal radiation heat exchange between the DHE and the guard vessel wall, 
which depends mainly on the temperature difference between these two 
components. 

Because the rise/down pipes and the reactor downcomer (adiabatic sections) 
are thermally insulated, there is no heat loss in the radial direction from these 
components and, therefore, U is equal to zero. The only heat exchange in these 
components is through the convection and conduction in the axial direction. 

To calculate the spatial distribution of the local coolant bulk temperature in 
the DHRL, each component in the loop is discretized into small spatial segments. 
Different components are discretized into different size segments. Non-uniform 
meshes are used for each component to keep the mesh size approximately the 
same at the junctions between components. In each segment, Eq. (22) is solved 
using a modified upwind finite difference method, which includes the 
contributions of both convection and conduction. 12 More detail will be discussed 
in the following sub-section. 

4.2.2 Numerical Method for Solving the Energy Balance 
Equation of the DHRL 

The DHRL is descritized into non-equal mesh size to accommodate the 
sudden change at the junctions between components. As shown in Fig. (14), the 
energy balance for the ith segment (control volume) can be expressed as: 



The left hand side of Eq. (26) is the specific heat change of the control 
volumn. In the right hand side, the first and the second terms are the net heat 
flux to the control volumn through axial convection and conduction and radial 
conduction respectively. Using the ’Power Law' scheme recommended by 
Patanka, 12 including the effects of convection and conduction, the heat flux 
density in the left and right boundaries of the ith segment, Ji and Ji+i, can be 

written as: 

J i = mCi x .-i + D^FtPj) fit 1 - it) 

A' s • (27a) 



(27b) 


where 

pi _ mCj 

6 A' s Dct . (27c) 

Equation (27) is applied throughout the DHRL, except at the entrance of 
the core, z=0, where Ji = Ji+i. 
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In these equations, the coefficiants D e r* and D e r i+1 are the effect thermal 
conductances in the left and right boundaries of the control volumn respectively. 
F(x) is a function considering dominate mode between convective and conductive 
heat transfer; 12 they are given as: 


D* — 2 


< 

+ 

i 

3 


kjr 1 kb , 

(27c) 

R 

it 

K> 




kb kb +1 , 

(27d) 


F(x) = Max (0,(1 - 0.1 |x| ) 5 ) + Max(0, -x). 


(27e) 


The function Max (a,b) is equal to the maximum value of a and b. 
Substituting Eqs. (27a) and (27b) into (26) and rearrange the results yields: 


dTt 

dt 


= KtuTj 1 +K U T{, +K u+1 n +1 +Qj 


i+1 


(28) 


Where N is the total number of control volumns in the DHRL. 


ir _ ihQ + De r AsF (Pi) 
JHi-l : 

MiCj, 


(28a) 


Ki,i = - 


DUsF (Pj) + mCj+i + D^ATF (Pj +1 ) + hjAj 

m4 


(28b) 


y. - DS‘A‘ +1 F (P* +1 ) 
1 , 1+1 : 

MiCj, 


(28c) 
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Qi = 

MiQ . (28d) 

Equation (22) was discretized into a set of algebra equations using a finite 
difference method. To solve these equations, Eq. (28) was in a matrix format as 
follows: 

{Tb} = [K] {T b > + {Q} (29) 


where: 




K U K,,2 
K2,l K2,2 

Kn,1 Kn,2 



The operator, [K] is a tri-diagonal matrix, except that Ki tN and Kn.i are 
not zero. 


To calculate the transient coolant bulk temperature in the DHRL, Tb, 
Trapezoidal time integrator was applied: 
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(30) 


(T b r 1 = (Tb] n + S [a (f b ) n+1 + (1- a) (f b } n . 

In Eq. (30) (t b P and |T b )" are obtained by using Eq. (29) in nth and 
(n+l)th time steps as: 

{f b | n = [Kf |T b )" +{Q} n > (31) 

(f b ) n+1 = [Kf* 1 |T b ) n *‘ +(Q) n *‘ , (32) 

combining Eqs. (30), (31) and (32), yields the following general form: 

( [I] - sa[Kf*‘) (T.,)"* 1 =( [I) + s(l-a) [Kf) [T b f 

+ s[a[Qf*' +(l-a)[Qj“] (33) 

Because thermal conductance matrix [K] is a function of coolant bulk 
temperature, Eq. (33) is a nonlinear equation. It was solved using a iterating 
method with fully implicit scheme ( a = 1.0 ). 


4.2.3 Hydraulics Model 


The coolant mass flow rate in the DHRL, by natural circulation after 
shutdown of the nuclear reactor, is determined from the overall momentum 
balance equation in the DHRL as follows [see Appendix B]: 



pg m A(z)dz - 



A(z)dP)oss 


(34) 


The integrations in Eq. (34) are carried out over the entire DHRL. The 
pressure loss term on the right hand side of Eq. (34) includes losses due to 
friction, pipe expansion and contraction, and acceleration losses of the coolant in 
the DHRL. 
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Fig. 15 Reactor Core Flow Channel and Wire Wrapped Fuel Pin. 
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In the reactor core region, all the flow channels are assumed to be identical 
as shown in Fig. (15). Therefore, the difference in flow rate between the center 
channels and the edge channels is negligible. The pressure losses in the core is 
calculated using the CRT (Chui-Rohsenow-Todreas) model for the triangular 
pitched channel with wire wrapped fuel rods: 1 3 

APffi = f s [ 1 +C = - AAP 2 

\2D e A^p c / (_ AHc( tcP 2 + H? )J , (35) 

where the friction factor, f s , is based on the pressure drop in the flow channel 
without a wire wrap, and 


D _ 2f}P 2 - K (Df 2 - s 2 ) 


K (Df + s ) 

(36) 

A c = SS. ( Df +s ), and A = £ P 2 - 7®? 
6 4 4- 

(37) 

In the pipes, reactor downcomer, and the DHE duct, 
calculated using Blasius relation for smooth walls: 

the friction losses are 

APE = X * 2 „ 

r 2 pw , 

(38) 

where the friction coefficient is given as: 


f s = (64/Re), for Re <2,100, 

(39a) 

f s = -0.001476 + 1.522 x 10 5 Re, for 2,100 < Re < 3,000, 

(39b) 

f s = 0.0058 + 0.05 / Re 0 - 32 , for 3,000 < Re < 3xl0 6 . 

(39c) 


The curve of f s vs Re is shown in Fig. (16). In this figure the boundaries of 
each flow regions are illustrated. 
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Eq. (29a) 
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Friction Factor as a Function of Reynolds Number. 




Equations (22) and (34) are solved iteratively to determine the coolant 
mass flow rate and the spatial distribution of the local coolant bulk temperature in 
the DHRL as functions of time after reactor shutdown. The numerical algorithm 
is explained in the block diagram of Fig. (17). 



A Block Diagram for the Algorithm of the DHRL Model. 

40 


Fig. 17 





4.3 DHE/Guard Vessel Heat Pipes Heat Rejection Model 


As indicated earlier, the DHE is radiatively coupled to the inside surface of 
the guard vessel. The heat removal from the DHE, and hence the coolability of 
the SP-100 reactor after shutdown, depends on both the energy transport 
capability of the guard vessel heat pipes to the decay heat radiators and the heat 
rejection capacity of the latter. At any time step in the transient numerical 
solution, the temperature of the guard vessel heat pipes, Thp, is determined from 
a quasi steady-state heat balance between the DHE and the guard vessel as: 



(T b (z,t) - T hp ) dz= Q re j, 


(40a) 


where, 

Qrej = h r A r ( Thp - T a ), or Psonic(Thp) , (40b) 

However, the temperature of the guard vessel heat pipes is determined by the 
dominant mode of heat transport, radiator heat rejection, Qrej = hr Ar (Thp ■ 
T a ), or the sonic limit of guard vessel heat pipes, Q re j = Psonic- 7,9 

At the sonic limit, the vapor flow at the exit of the evaporator region of the 
heat pipe is choked, where the flow velocity is equal to the sonic speed of the 
vapor at the evaporator temperature; T hp- The maximum power throughput due 

to the sonic limit is given by: 23 

Psonic = 0.474 Ahphfg V Pv pv , (40c) 

where hf g , p v and p v for the heat pipe working fluid are evaluated at Thp. 

Initially, following the reactor shutdown, the guard vessel heat pipes 
temperature is such that heat rejection by thermal radiation from the decay heat 
radiators to the ambient on the lunar surface is lower than the energy transport 
by the vessel heat pipes at the sonic limit. In this case, the decay heat removal 
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from the reactor core is restricted by the surface area of the decay heat radiators 
(see Fig. (18)). However, as the decay heat power decreases with time, both the 
guard vessel wall temperature and the sonic limit of the vessel heat pipes will 
decrease. Eventually, as the sonic limit of the heat pipe drops below that of the 
decay heat radiator, the heat removal from the DHRL will be restricted by the 
sonic limit of the heat pipes. That will slow down the heat removal from the SP- 
100 reactor core (see Fig. (18)). 


Although the heat pipe sonic limit would slow down the decay heat removal 
from the reactor core, it could prolong the time for the lithium coolant to cool 
down to its freezing temperature. 14 Lithium freezing in the primary loop is not 
desirable because the formation of voids during freezing could induce hot-spots 
in the reactor core during a subsequent startup.15-17 Therefore, guard vessel heat 
pipes sonic limit can be used to prevent the lithium freezing with a sacrifice on 
the decay heat removal capability of the DHRL system. The designers must make 
some compromise in these two aspects. 

For different heat pipe working fluids (water, potassium and sodium) and 
heat pipe diameter, the heat pipe code HTPIPE20i s used to calculate the sonic 
limits for the guard vessel heat pipes at the evaporator temperature and the 
results are presented in Fig. (18). The heat rejected from the decay heat radiator 
with a different surface area is calculated as a function of the radiator 
temperature and the results are also present in Fig. (18). We assumed the 
temperatures of the guard vessel heat pipes and of the decay heat rejection 
radiator are equal. As delineated in this figure, for a decay heat radiator surface 
area of 25 m2, the potassium heat pipes of 2.54 cm diameter in the guard vessel 
heat pipes reach their sonic limit at an evaporator temperature of 111 K. This 
transition temperature, from a radiator limited heat removal to a heat pipe 
limited heat removal, increases as the radiator area increases. Fig. (18) also 
demonstrates that water heat pipes are unsuitable for cooling the guard vessel 
walls because the wall temperature would initially exceed the critical temperature 
of water (~ 550 K) for a medium decay heat radiator (such as 25 cm2, see t h e 
discussion in Section 5). Therefore, potassium heat pipes are a better choice for 
cooling the guard vessel wall. The effect of the transition from a radiator limited 
heat rejection to heat pipes limited heat rejection on the coolability of the DHRL 
is included in the decay heat removal model. In our analysis, potassium heat 
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its of Diameter and Working Fluid of the Guard Vessel Heat 
i on Decay Heat Rejection. 








pipes are used to transport the heat from the guard vessel wall to the decay heat 
radiator on the lunar surface. 

In this section, the transient overall momentum and energy balance 
equations of the DHRL described, together with the transient heat conduction 
equation in the fuel pin are solved interactively using a fully implicit time 
integrator. The coupled decay heat removal model calculates the spatial 
distribution of the local bulk temperature of Li coolant in the DHRL as well as 
the axial and radial temperature distributions in the fuel pin, as functions of time 
after reactor shutdown. In the following section, the results of the analyses 
investigating the coolability of the SP-100 reactor core after shutdown by natural 
circulation on the lunar surface are presented for the base case design and 
operation parameters listed in Table 1 . The parametric analyses of the effects of 
the design parameters of the DHRL on the decay heat removal capability of the 
system are also performed and the results are included in the next section. 
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5. RESULTS AND DISCUSSION 


This section includes three parts. The first part presents the results of some 
transient analyses of the passive decay heat removal capability of the DHRL of the 
SP-100, 550- kW e power system for a lunar outpost following a LOFA. In the 
second part, the results of the parametric analyses for some design parameters of 
the DHRL are presented and discussed. In the last part, the differences and 
similarities of the thermal behavior of the system after shutdown at two gravity 
constants (lg and 1/6 g) are discussed. These results would be useful in relating 
earth-based preflight test results to the actual system performance on the lunar 
surface. 

5.1 Base Case Result 

The base case design and operation parameters for the SP-100, 500-kWe 
power system for a lunar outpost, listed in Table 1., are based on the information 
provided by Rockwell International Corp. 22 The envelop dimensions of the 
reactor and the shield are shown in Fig. (19). The nominal thermal power of the 
SP-100 reactor is about 2.3 MW thermal, at which the coolant temperature at the 
exit of the reactor core is ~ 1355 K. The reactor core consists of approximately 
1296 fuel pins in a triangular arrangement. 

The initial condition for our calculation is based on the following 
hypothesis: At a time when the SP-100 reactor is operating at its nominal power 
level, the reactor scrams, and at the same time, the lithium coolant mass flow rate 
drops to zero (LOFA). Therefore, the vanturi allows natural circulation of 
lithium coolant through the DHE initiated by the density difference between the 
rise and down pipes. To be practical, the coolant mass flow rate may not drop to 
zero at the same time of reactor scram. These initial conditions should be 
considered as the worst case in a LOFA. 
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5.1.1 Spatial Coolant Temperature Distribution 

Figure (20) presents the calculated spatial distributions of the local bulk 
temperatures of Li coolant in the DHRL at different times after reactor 
shutdown. The steep temperature gradients at the interfaces between the reactor 
core, both the adiabatic rise and down pipes, and between these pipes and the 
DHE, clearly demonstrate that the effect of axial heat conduction in the Li coolant 
flow is important for the heat transport in the DHRL by natural circulation. As 
indicated in Fig. (20), immediately after reactor shutdown, the maximum coolant 
temperature in the DHRL (~ 1355 K) occurs at the exit of the reactor core. This 
is also the case during full power operation, where forced convection is the 
dominant mode of heat transfer. However, since the coolant flow rate by natural 
circulation is much lower than that at full power operation (see Fig. (22) and 
Table 1.), the contribution of axial heat conduction in the liquid metal coolant, 
immediately after shutdown, becomes important. The axial heat conduction in Li 
coolant, together with convective heat transfer, causes the maximum coolant 
temperature in the DHRL to increase (up to 100 seconds after shutdown) and shift 
with time away from the exit of the SP-100 reactor core toward the DHE. 

As shown in Fig. (20), it takes about 200 seconds for the hottest part of the 
Li coolant (with the maximum coolant temperature) to travel from the exit of the 
reactor core to the entrance of the DHE. In addition to being shifted toward the 
DHE, the temperature of the hottest part of Li in the DHRL decreases with time 
after shutdown from a peak value of 1530 K, happened at about 50 seconds after 
shutdown, to about 1250 K after 300 seconds. 

The results in Fig. (20) also show that at times greater than 200 seconds 
after shutdown, the maximum coolant temperature in the DHRL occurs at the 
entrance of the DHE, while the lowest coolant temperature always occurs at the 
exit of the DHE. This is because the DHE acts as a heat sink in the DHRL after 
the reactor shutdown. The lowest coolant temperature also decreases from 
approximately 1275 K, immediately after shutdown, to about 900 K after 300 
seconds. Despite the rapid decrease in decay power in the reactor core with time 
after shutdown, the average coolant temperature in the DHE decreases much 
slower. Such a slow decrease in the coolant temperature in the DHRL after 


47 



300s 

200s 



(») 3UniVH3dW31 1NV1000 


48 


Fig. 20 Spatial Distribution of the Lithium Coolant in the DHRL After 
Reactor Shutdown on the Lunar Surface. 



reactor shutdown is because of the large heat capacity of the lithium coolant ( > 
4.5 kJ/kg K) and a very low coolant mass flow rate. Although the maximum 
coolant temperature in the DHRL does not drop below its value at full power 
operation (-1355 K, see Table 1.) before 50 seconds after shutdown, the coolant 
temperature in the reactor core insignificantly increases above its value during 
full power operation, instead it drops very rapidly after shutdown. 

It is worth noting that during the first 100 seconds (even more) after the 
reactor shutdown, the coolant temperatures in the rise pipe is higher than its 
value during nominal full operation. Sometimes (like - 50 seconds after 
shutdown), its maximum value is as much as 175 K higher than its nominal value. 
This temperature increase might contribute to some undesired effects (such as 
boiling in the lithium coolant or creep in the pipe walls) in the system. 

Figure (21) shows the axial variation in the centerline temperature of a fuel 
pin at the center of the reactor core after shutdown. As this figure indicates, 
after shutdown, the maximum fuel temperature always occurs at the exit of the 
reactor core and decreases to a relatively lower temperature (< 1200 K) after 300 
seconds following reactor shutdown. 

5.1.2 Transient Maximum Coolant and Fuel Pin Temperatures 

Figure (22) compares the calculated maximum coolant temperatures in the 
DHRL with those of the fuel and cladding in the SP-100 reactor core after 
shutdown. Also delineated in Fig. (22) are the calculated time dependent 
temperatures of the guard vessel (or its heat pipes) and the coolant mass flow rate 
by natural circulation in the DHRL. As this figure shows, the maximum fuel 
temperature drops very rapidly from its nominal value of about 1950 K at full 
power operation (see Table L), to about 1500 K within the first 5 seconds after 
shutdown. At such time, the very low coolant flow rate by natural circulation 
results in the removal rate of the decay heat from the fuel pins being lower than 
its rate of generation, causing the maximum fuel temperature, at the centerline of 
the fuel pin to rise reaching a maximum of -1600 K at 30 seconds after 
shutdown. Beyond this point, the decay heat generation rate drops below the rate 
of removal by natural circulation, causing the maximum fuel temperature to 
decrease with time, reaching as low as 1000 K at 1000 seconds after shutdown. 
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Fig. 21 Axial Distribution of the Fuel Pin Centerline Temperature as 
Function of Time After Shutdown. 
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Fig. 22 Calculated Maximum Coolant, Fuel, and Cladding Temperatures and 
Coolant Mass Flow Rate in the DHRL as Functions of Time after 
Shutdown. 




Results in Fig. (22) indicate that, while the maximum temperature of the 
coolant in the DHRL peaks shortly after shutdown (~ 30 seconds), the coolant 
mass flow rate by natural circulation takes up to 120 seconds to reach a maximum 
value of -0.14 kg/s. This increase in the coolant mass flow rate occurs when the 
hot coolant from the reactor core travels through the adiabatic rise pipe(see Fig. 
(7)). The induced density difference between the hot liquid in the rise pipe and 
the cold liquid in the down pipe causes the coolant mass flow rate in the DHRL to 
increase with time. It reaches a maximum value of 0.14 kg/s as the hot coolant 
enters the DHE duct. When the hot coolant exits the rise pipe and travels through 
the DHE duct, its temperature decreases and hence the buoyant force for natural 
circulation, causes the coolant mass flow rate to decrease. Such a decrease in the 
coolant mass flow rate continues until the coolant from the DHE begins to travel 
through the down pipe toward the reactor core. This occurs at about 400 seconds 
after shutdown, at which time the coolant mass flow rate reaches a minimum 
value of - 0.075 kg/s. When the coolant from the DHE fills the down pipe, the 
density difference between the hot coolant in the rise and the cold coolant in the 
down pipe increases, causing the coolant mass flow rate to increase again, it 
reaches a second, but much lower, peak value of ~ 0.082 kg/s, at approximately 
520 seconds after shutdown. This response of the coolant mass flow rate 
continues each time when the cold coolant from the DHE travels through the 
down pipe. However, the peak mass flow rate decreases progressively with time 
after shutdown, commensurate with the decrease in decay power. Such a 
variation in the coolant mass flow rate in the DHRL is responsible for the 
variations shown in the coolant, fuel, and cladding temperatures delineated in Fig. 
( 22 ). 


While the responses of the maximum cladding and coolant temperatures are 
similar to that of the maximum fuel temperature, the temperature of the guard 
vessel wall generally decreases very slowly after shutdown (see Fig. (22)). Such 
slow decrease in the guard vessel heat pipe temperature, particularly at times in 
excess of -120 seconds after shutdown, is because the decay heat removal from 
the DHRL is constrained by the sonic limit of the guard vessel heat pipes. 
Conversely, during the first 120 seconds after shutdown the decay heat removal 
from the DHRL is constrained by the decay heat radiator (surface area and 
temperature) on the lunar surface, causing the guard vessel wall temperature to 
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decrease at a relatively higher rate. This thermal behavior will be discussed in the 
second part of this section. 


The above-mentioned results suggest that natural circulation of lithium 
coolant in the DHRL of the SP-100, 550-kWe power system for a lunar outpost 
would maintain the SP-100 reactor core safely coolable after shutdown. 
However, the lithium coolant in the adiabatic rise pipe, directly downstream of 
the reactor core, could be overheated by as much as 175 K, above its nominal 
operation value, for approximately 200 seconds after shutdown. The effect of 
such temperature overheating on the structure strength of the rise pipe walls 
should be assessed before a final conclusion can be drawn regarding the 
coolability of the DHRL by natural circulation of lithium coolant. In the second 
part of this section, the effects of changing some design parameters of the DHRL 
on improving the coolability of the DHRL are being investigated. 

5.2 Parametric Analyses of the Coolability of the DHRL 

In this sub-section, the decay heat removal model is used to examine the 
effects of the following parameters on the decay heat removal capability of the 
system: (a) the height of the DHE duct; (b) the elevation of DHE, (c) the diameter 
of DHRL pipes; (d) the surface area of the decay heat rejection radiator on the 
lunar surface; and (e) the diameter of the guard vessel heat pipes. In the previous 
sub-section, it was shown that the maximum coolant temperature of the DHRL 
could be as much as 175 K above its nominal value with the decay heat removed 
by natural circulation after reactor shutdown. Other components' temperatures 
are found to be below their nominal values during full power operation. 
Therefore, the maximum coolant temperature after shutdown will be taken as the 
criterion for assessing the coolability of the system. 

5.2.1 Effect of DHE Dimensions 

Figure (23) presents the results on the effects of changing the height of the 
DHE duct, with the same aspect ratio (see Table L), and the surface area of the 
decay heat radiator on the coolability of the DHRL by natural circulation after 
shutdown. In the base case, the height of the DHE rectangular duct is 10 cm (see 
Table 1.), for which the maximum coolant temperature in the DHRL exceeds 
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1550 K after reactor shutdown. However, in order for this temperature to be 
brought closer to its value during nominal operation of the power system (1355 
K), the height of the DHE is increased. In addition to increasing the surface area 
of DHE, such increase in the duct height reduces friction losses, resulting in a 
higher coolant flow rate and lower coolant temperature. Conversely, increasing 
the height of the DHE duct beyond 20 cm, increases friction pressure losses, 
hence reducing the mass flow rate. This increase in friction losses is because the 
coolant flow will be in the transitional regime, where 2100< Re < 3000 (see Eq. 
(39b)). Results shown in Fig.(23) also indicate that increasing the height of the 
DHE duct beyond 15 cm, the resulting decrease in the maximum coolant 
temperature is very small, hence that the change in the DHE duct height, when it 
is beyond 15 cm, will insignificantly affect the coolability of the DHRL. At such 
a duct height, the maximum coolant temperature in the DHRL after reactor 
shutdown is about 1480 K, which is still about 125 K higher than that at full 
poweroperation. 
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Rg. 23 Effect of Dimensions of the DHE Duct on the Coolability of 
the SP-100 Reactor after Shutdown. 




Results show that even if the DHE duct height is doubled, the maximum coolant 
temperature after shutdown will still be more than 100 K higher than its nominal, 
full power operation value. 

Figure (23) also indicates that increasing the surface area of the decay heat 
radiator on the lunar surface has no effect on the maximum coolant temperature, 
but reduces the radiator temperature, and, to a less extent, the coolant flow rate 
in the DHRL. 

5.2.2 Effect of Elevation of DHE and Diameter of DHRL Pipes 

Figure (24) presents the results on the effects of changing the elevation of 
the DHE and/or the diameter of the rise and down pipes in the DHRL on the 
coolability of the SP-100 by natural circulation after shutdown. As this figure 
shows, increasing the elevation of the DHE or increasing the diameter of the 
pipes improves the decay heat removal capability of the DHRL. However, 
increasing the diameter of the rise and down pipes beyond 8.5 cm insignificantly 
affects the maximum coolant temperature in the DHRL. The results in Fig. (24) 
demonstrate that increasing the elevation of the DHE increases the circulation rate 
of the coolant, resulting in a relatively lower maximum coolant temperature in 
the DHRL after reactor shutdown. For example, for a pipe diameter of 8.5 cm, 
increasing the elevation of the DHE from 1 .96 m to 2.94 m increases the coolant 
flow rate by about 24 % (from 0.125 to 0.155 kg/s), but lowers the maximum 
coolant temperature in the DHRL by only 36 K (from 1562 K to 1526 K). 
Therefore, increasing the elevation of the DHE to improve the decay heat 
removal capability of the DHRL is not recommended because of the added system 
mass and excavation requirement. Also, increasing the diameter of the rise/down 
pipes is not recommended since it will affect the reactor design and the secondary 
loop. Instead, increasing the the height of the DHE duct, as discussed in the 
previous section, could better enhance the coolability of the SP-100 reactor after 
shutdown by natural circulation. 
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5.2.3 Effects of Radiator Area and Guard Vessel Heat Pipe 
Diameter 

Figures (25a) and (25b) show the effect of the surface area of the decay 
heat rejection radiator and the diameter of the guard vessel heat pipes on the 
guard vessel wall temperature and the maximum coolant temperature after 
reactor shutdown. Fig. (25a) shows the base case guard vessel heat pipe diameter 
of 2.54 cm, while Fig. (25b) is for a larger heat pipe diameter of 3.81 cm. As 
delineated in these figures, within the first 20 seconds after shut down, the 
maximum coolant temperature is independent on the radiator area and/or the 
diameter of the guard vessel heat pipes. However, at later times up to 200 seconds 
after shutdown, increasing the radiator area from 12.5 to 25 m 2 , slightly lowers 
the maximum coolant temperature. Further increase in the decay heat rejection 
radiator to 50 m 2 does not affect the maximum coolant temperature, hence the 
coolability of the system. It can also be seen in these figures that at times longer 
than 200 seconds after shutdown, both the radiator area and the diameter of the 
guard vessel heat pipes insignificantly affects the coolant flow rate or the 
maximum coolant temperature. 

Figures (25a) and (25b) also show that increasing the guard vessel heat pipe 
diameter does not affect the maximum coolant temperature, but reduces the guard 
vessel temperature and increases the time after shutdown, beyond which decay 
heat rejection is restricted by the sonic limit of the heat pipes (see Section 4.3). 
For a heat pipe diameter of 2.54 cm, the sonic limit is reached at 20 seconds and 
110 seconds after shutdown for a radiator area of 25 m 2 and 12.5 m 2 , 
respectively (see Fig. (25a)). For the largest radiator area of 50 m 2 , the radiator 
temperature (equal to the guard vessel temperature) immediately after shutdown 
is already below the sonic limit of the guard vessel heat pipes (see Figs. (25a) and 
(25b)). When the heat pipe diameter in the guard vessel is increased from 2.54 
cm to 3.81 cm, the time to reach the sonic limit increased to 45 seconds and 190 
seconds for a radiator area of 25 m 2 and 12.5 m 2 , respectively. Therefore, using 
a larger guard vessel heat pipe can defer the time for reaching the sonic limit for 
the same radiator area. 
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Fig. 24 Effect of Elevation of the DHE and Diameter of Rise and Down 

Pipes of the DHRL on the Coolability of the SP-100 after Shutdown. 





5.3 Effect of Gravity on the Decay Heat Removal 

To identify the differences and similarities between test results at earth 
gravity and the actual system performance on the lunar surface, the effect of 
increasing the gravity constant from 1/6 g on the lunar surface to lg on the 
thermal behavior of the DHRL was investigated. The results presented in Figs. 
(26)-(28) indicate that at earth gravity, the decay heat removal capability of the 
power system is significantly higher. As shown in Fig. (26), unlike on the lunar 
surface (see Fig. (20)), the maximum coolant temperature does not exceed its 
nominal value (at t=0). That means that in an earth-based preflight test, there 
will be no coolant overheat happening during natural circulation cooling after the 
reactor shutdown, even though this may happen on the lunar surface. 

Figure (27) shows that the maximum coolant temperature in the DHRL at 
lg, immediately after shutdown, is approximately 1380 K, which is about 140 K 
lower than its predicted value on the lunar surface. Moreover, the maximum fuel 
and cladding temperature could be as many as 50 K to 100 K lower than their 
values on the lunar surface (see Fig. (22)). Also, the maximum fuel temperature 
in the SP-100 reactor after shutdown would be lower, but the coolant flow rate 
would be higher than their values on the lunar surface. As Fig. (28) shows, the 
coolant flow rate by natural circulation in the DHRL during earth-based testing is 
more than twice that on the lunar surface. Furthermore, the frequency of the 
coolant flow oscillation on earth is higher than on the lunar surface, for the 
coolant mass flow rate on 1 g is higher. 

These comparisons show that there is a large difference between the 
thermal behavior of the system after the reactor shutdown at 1 g and on the lunar 
surface. Recognizing these differences, an earth-based preflight test results could 
be used to predict the actual system performance on the lunar surface. 
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Calculated guard vessel and maximum coolant temperatures for 
diameter d=2.54 cm guard vessel HP and radiator area. 
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Fig. 25b Calculated guard vessel and maximum coolant temperatures for 
diameter d=3.81 cm guard vessel HP and radiator area. 
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6. SUMMARY AND CONCLUSIONS 

The decay heat removal from an SP-100 reactor of a 550 kWe power 
system for a lunar outpost by natural circulation of lithium coolant was 
investigated. A transient model, which simulates the DHRL of the power system 
and incorporates axial conduction in the liquid metal coolant, has been developed 
and used to assess the system’s decay heat removal capability. The effects of the 
surface area of the decay heat rejection radiator, dimensions of the DHE flow 
duct, elevation of the DHE, guard vessel heat pipes diameter, diameter of DHRL 
pipes, and gravity on the decay heat removal capability were examined. 

Results showed that natural circulation of the lithium coolant in the DHRL 
would maintain the SP-100 reactor core safely coolable after shutdown. 
However, the lithium coolant in the adiabatic rise pipe, directly downstream of 
the reactor core, could overheat by as much as 175 K above its nominal operation 
value of 1355 K at approximately 200 seconds after shutdown. Such a coolant 
temperature overheat, which lasted for a very short time, can be reduced by as 
much as 50 K by increasing the height of the DHE duct up to 15 cm. A further 
increase in the DHE duct height, would have little effect on the decay heat 
removal from the SP-100 reactor core after shutdown. Increasing the elevation 
of the DHE slightly improves the decay heat removal capability of the DHRL, but 
does not justify the additional system mass and excavation on the lunar surface. 
Therefore, increasing the elevation of the DHE is not recommended since it does 
not significantly improve the system’s decay heat removal coolability. It would 
strongly impact the mass and size of the system, and hence the excavation and 

launch costs. 

Increasing the area of the decay heat radiator insignificantly affects the 
coolability of the DHRL; a radiator area in the order of 10-15 m 2 , will be 
sufficient to maintain long term coolability of the power system by natural 
circulation. Analysis demonstrated that increasing the diameter of the heat pipes 
in the guard vessel does not influence the coolability of DHRL, it lowers the 
guard vessel wall temperature. Increasing the diameter of DHRL pipes up to a 
value of 8.5 cm insignificantly affects the coolability of the system, and therefore 

it is not recommended. 
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TEMPERATURE (K) 


Results also showed that the thermal performance of the power system 
after reactor shutdown in an earth-based preflight test would be significantly 
different than that on the lunar surface. The maximum coolant temperature in the 
DHRL and the maximum fuel temperature in the reactor core could be as muc 
as 140 K and 50 - 100 K lower than their values on the lunar surface. These 
results would be useful for relating data earth-based tests results to the decay heat 
removal capability of the actual power system on the lunar surface. 



Fig. 27 Effect of Gravity on Calculated Maximum Coolant and Fuel 
Temperatures after Reactor Shutdown. 
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Effect of Gravity on the Coolant Mass Flow Rate in the DHRL after 
Reactor Shutdown. 




Table 1. Base Case Parameters for an SP-100,500-kWe Power System 
for a Lunar Outpost. 


Parameter 


Value 


1. EXTERNAL PIPES AND DHE 

Pipe height (m) 

Pipe inner diameter (m) 

Aspect ratio of DHE duct 
Height of DHE duct (m) 

Diameter of DHE circle (m) 

Pipes wall thickness (m) 

DHE Na heat pipe length (m) 

Number of DHE Na heat pipes 
Diameter of guard vessel heat pipes (m) 
Pipe wall material 
DHE pipe surface emissivity 

2. CHARD VESSEL AND DECAY 

Guard vessel height (m) 

Guard vessel inner diameter (m) 

Number of guard vessel heat pipes 
DHE/guard vessel radiation view factor 
Emissivity of DHE surface 
Emissivity of guard vessel wall 
Radiator surface area (m 2 ) 


2.45 

0.107 

0.26 

10.0 

2.07 

0.005 

0.3 

100 

2.75 

Nb-1% Zr 

0.8 

HEAT RADIATOR 

3.53 

2.30 

45 

1.0 

0.8 

0.8 

25 


V REACTOR CORE 

Effective core height (m) 

Height of core vessel (m) 

Reactor core effective radius (m) 
Effective inner radius of core vessel (m) 
Total no. of fuel elements 
Fuel rod radius (m) 

Fuel rod lattice 
Pitch-to-diameter ratio (P/D) 
Fuel-cladding gap size (mm) 

Cladding material 
Fuel material 

Diameter of wire wrap (mm) 

Coolant type 


0.3175 

0.4258 

0.1764 

0.2546 

1296 

0.0037 

triangular 

1.07 

0.13 

Nb-1% Zr 
UN 
0.0548 
Lithium 


4. OPERATION PARAMETERS BEFORE SHUTPOm 

Reactor thermal power (MW) 

Core inlet coolant temperature (K) 

Core exit coolant temperature (K) 

Core coolant mass flow rate (kg/s) 


2.3 

1260 

1355 

14.1 
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APPENDIX A 


Estimation of the DHE/Guard Vessel View Factor 

The view factor between DHE and guard vessel was calculated according 
the geometries and dimensions of these two components. After simplification, the 
configuration of DHE/guard vessel is shown in Fig. (A-l). The DHE is a circular 
ring identified by 'L' in Fig. (A-l). The guard vessel is a hollow right column 
identified by ' 3 ' in Fig. (A-l). According to reference 6, the view factor 
between the DHE and Guard Vessel, F1.3, is calculated by. 

Fi.3 = £ [ L + bF(b) + cF(c) - (L + b) F(L+b) - (L + c) F(L+c)] , (A-l) 

where: 


(A- 2 ) 



b, c, L, r and R are defined in Fig. (A-l). 

In our calculation, the DHE is discretized into some small segments. In the 
jth segment as shown in Fig. (A-l), the radiation view factor with the guard 
vessel , Fj-3, should be equal to the view factor between the DHE and guard 

vessel, F1-3: i.e. 

F 1-3 = Fj-3. 
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Pig. ( A - 1 ) Configuration of DHE/guard vessel. 
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Vessel/ 


By the definition of view factor: 
A1F1-3 = A3 F3-1, 


(A- 3 ) 


where A, and A 3 are the surface area the DHE and guard vessel respectively. For 
jth segment, also by definition: 


A3F3-J - Aj Fj. 3> 


(A- 4 ) 


and 


F3-1 = X Fj-i 
j 

Substituting Eqs. (A- 4 ) and (A- 5 ) into Eq. (A- 3 ) yields: 

A1F1-3 = £ AjFj.3 
i 

Because 

a, = X a j 

j 

and Fj. 3 is a constant for a circular DHE, 

F1-3 = Fj -3 • 


(A- 5 ) 


(A- 6 ) 


(A- 6 ) 
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APPENDIX B 


Momentum Balance Equation for DHRL 

For the flow channel shown in Fig. (B-l), a control volumn, R, is selected. 
Assume that between time, t, and time, t+At, this control volumn changes from 
R= Ri + R2 to R'= R.2 + ^3, causing the fluid momentum to ^ 

M r = M r1 +Mr2 to Mr. =M r2 +M b3 , where the vector M E ' and M m represent 
the fluid momentum in the control volumn R2 at t+At_ and t, respectively. 
Therefore, during the interval At, the rate of change of the fluid momentum can 
be expressed as: 

AM r Mr. -M r _ Mr2' - Mr 3 - M R1 
At ~ At At At , (B-l) 

The first term on the RHS of Eq. (B-l) is the momentum change in R2; as 
At— »0,S'— »Sand R'— » R(see Fig. (B-l)), this term can be expressed as: 


lim 

At— *0 


Mr2’ " Mr2 


At 



(B-2) 


where pR is the fluid density, V is its velocity vector in the control volumn R. 

Similarly, the second term on the RHS of Eq. (B-l), which is the 
momentum change between volumn Ri and R3, can be expressed as: 


lim 
At — *0 




At 


-1 


V s dm 


(B-3) 
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where V s is the fluid velocity vector at the surface S. The mass flow rate in the 
control volumn, R, can be expressed as. 

dm=p s v s dS <B ' 4) 

Substituting Eq. (B-4) in to Eq. (B-3) and rearranging the result, yields: 


lim 

At-*) 




At 


= ((PsVs)v s dS 


(B-5) 


From the Newton law, the total force acting on the control volumn, F, is 
equal to the rate of change in momentum within the control volumn: 

F=jf PR VdR + / M V s dS (B-6) 

Ot*R /s 

For one-dimensional slug flow (i.e. the V s is uniform at the inlet (S i) and 
outlet (S2) surfaces of the control volumn), the second term of Eq. (B-6) can be 
simplified as: 

J s (ps ^s)'V d S = J a '^ dm =m 2 V 2 -m 1 V 1 (B . 7) 

Substituting Eq. (B-7) intoto Eq. (B-6), the one-dimensional flow 
momentum equation is: 

p = i(MV) + m 2 V 2 -m 1 V 1 (B-8) 

dt 


72 


Through the use of flow continuity ( =m ) and substitute the mass of the 
control volume M = p AV, and the total force, F = ^ F i, Eq- ( B ' 8) can be 
presented as: 


AV — (pV)=X F i -m(V 2 -V 1 ) 
dt , 


(B-9) 


The DHRL, shown in Fig. (B-2), was modeled using Eq. (B-9) The 
whole loop was divided into N segments and the energy balance equation 
wassolving in each segment (see Section 4.2.2). In the id, segm^t except *ose 
located at the interface between different sections in the loop, Eq. (B-9) can be 

presented as: 


— — (pVA) = X F r m (V 2 -V,) 
A dt ; 


where A is the cross section area 
expressed as: 


(B-10) 

of the DHRL. The force balance in the loop is 


51 Fj = ± AV pg + ^ 2^2 ‘ Fi Aj - APfcfcA 

i 

Substituting Eq. (B-ll) into Eq. (B-10) and, AV = A AL, yields: 


(B-ll) 


AL l (m ) = ± ApgAL + P 2 A 2 - P, A , - m(V 2 - V ,) - AP tlic A 
dt 


(B-12) 


Because A is not the same in the different sections of the DHRL, Eq. (B-12) was 
integrated over the whole DHRL as follow. 


73 



A : Surface Area of the Segment 
AV : Volume of the Segment 
P : Pressure 
V : Coolant Speed 


Fig. (B-2). A Simplified Schematic of the DHRL. 
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| il(m) dz = | pgA(z)dz + 


dPA(z) - 


A(z) dPfnc + 


I. 


where L is the length of the DHRL. Since for the whole loop, 

[ d(PA(z)) =0 and J mdV(z) = 0 
Jh ’ 

then the overall momentum balance can be written as: 


A(z) dPfric 


| ^(m) dz = | pgA(z)dz - 
which is the same equation as Eq. (34) in Section 4.2.3. 


m dV(z) 

f 

(B-13) 


(B-14) 
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